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THE  STEERING  PROBLEM  FOR  THE  MIRROR-SCAN  TRACKING  SYSTEM 
PART  1 - THE  STOP-GO-STOP  PROCEDURE 


1.  INTRODUCTION 

The  Mirror  Scan  system  employs  a mirror  to  steer  the  beam  of  a tracking  radar  whose 
light  wei^t  (about  11  kg  (25  lb))  permits  it  to  be  slewed  much  more  rapidly  than  the 
antenna  of  a conventional  system,  so  that  the  system  can  be  used  to  track  several  targets 
simultaneously.  Since  tracking  accuracy  varies  inversely  with  update  time,  the  number  of 
targets  which  can  be  simultaneously  tracked  depends  on  the  efficiency  of  the  steering 
algorithm.  The  steering  problem  is  to  determine  the  radar-beam  motion  which  covers  a 
given  collection  of  targets  with  minimum  update  time,  subject  to  two  constraints: 

• A physical  constraint  on  the  magnitudes  of  the  torques  that  can  be  applied  to  the 
mirror; 


• A constraint  on  the  velocity  of  the  beam  motion  as  the  beam  passes  through  a 
target,  so  as  to  insure  a required  number  of  radar-pulse  hits  on  the  target. 


1.1  General  Synopsis 

Sections  2 through  4 are  devoted  to  a careful  mathematical  description  of  the  steering 
problem,  which  is  a difficult  problem  in  control  theory.  As  a first  approach  to  this  problem 
we  will  in  this  report  examine  the  stop-go-stop  (SGS)  paths,  that  is,  beam  motions  in 
which  the  beam  is  made  to  stop  at  each  target.  The  problem  of  finding  the  optimal  SGS 
p«th  can  be  reduced  to  a standard  problem  in  control  theory  with  a known  solution.  This 
reduction  of  the  problem  is  described  in  section  5.  In  section  7 we  wiU  show  how  calcula- 
tions of  update  times  of  optimal  SGS  paths  for  various  configurations  of  target  positions 
lead  to  the  conclusion  that  the  mirror-scan  system  can  track  at  least  seven  targets  in  a half- 
hemisphere with  an  update  time  of  less  than  1 second. 

In  section  6 we  will  discuss  the  aiming  problem , which  is  the  problem  of  determining 
the  best  system  orientation  for  a given  area  of  coverage.  One  mi^t  expect  that  the  solution 
of  this  problem  would  be  to  point  the  system  at  the  centroid  of  the  area  of  coverage,  and 
computer  calculations  suggest  that  this  is  in  fact  the  case.  However,  the  same  calculations 
also  reveal  that  the  system  orientation  has  only  a sli^t  effect  on  efficiency,  provided  that 
the  departure  of  the  pointing  direction  from  the  centroid  of  the  area  of  coverage  is  not  too 
“extreme.” 

Finally,  in  some  situations  the  SGS  procedure  is  nearly  optimal,  and  in  others  SGS  is 
clearly  inefficient.  In  subsequent  reports  in  this  series  we  hope  to  obtain  lower  bounds  on 
update  time  (as  a function  of  number  of  targets  and  area  of  coverage),  identify  the  precise 
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conditions  under  which  SGS  is  grossly  inefficient,  and  construct  better  tracking  algorithms 
for  this  case. 

1.2  Tedinical  Summary 

A smnmary  of  the  technical  results  of  this  report  is  as  follows : 

• The  target  position  ( and  the  corresponding  mirror  normal  position  p are  rdated  by 

{-a-2(a-p)p,  (2.1) 

sdiere  a is  a constant  vector  defining  the  orientation  of  the  system.  (The  aiming  vector  is 
-a.) 

• The  syetem  control  is  defined  by  specifying  the  applied  torque  vector  u as  a func- 
tion of  time. 

• The  Systran  control  u = u(I)  and  mirror  motion  p ^ p(0  are  related  by 

u = /pXp,  (3.1) 

where  / is  a certain  moment  of  inertia. 

• The  steering  problem  is  to  find  the  closed  curve  p = p(0  which  joins  given  points 
with  minimum  update  time,  subject  to  two  constraints: 

IpXpKa  (4.C1) 

and 

(4.C2) 

• Ihe  constraint  (4.C1)  can  be  expressed  as  a bound  on  the  second  covariant  deriva- 
tive 5^p/5l^  of  p,  which  is  the  tangential  component  of  p. 

• The  constraint  (4.C2)  is  closely  approximated  in  a more  convenient  form  by 

|pl<P/2.  (4.CA2) 

• In  the  stop-go-stop  (SGS)  procedure  the  constraint  (4.C2)  is  automatically 
satisfied,  and  p moves  in  great-circle  arcs,  stopping  at  each  “image”  of  a target  position. 

• The  determination  of  the  optimal  point-to-point  motion  of  p along  an  SGS  path 
reduces  to  a solved  problem  in  control  theory,  and  the  update  time  T of  an  SGS  path  is 
then  given  by 
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• There  are  (N  - l)!/2  SGS  paths,  where  N is  the  number  of  targets.  The  optimal  SGS 
path  is  the  one  which  minimizes  £ y/Zj  (and  therefore  7)  and  can  be  found  by  strai^tfor- 
ward  enumeration  of  all  cases,  provided  that  N is  not  too  large. 

• For  the  current  design  value  of  a it  can  be  predicted  from  the  calculated  results 
that  the  Mirror  Scan  system  can  simultaneously  track  at  least  seven  targets  in  a half- 
hemisphere with  an  update  time  of  less  than  1 second. 

• The  orientation  of  the  system  (defined  by  a)  has  only  a sli^t  effect  on  system 
efficiency,  provided  that  the  departure  of  the  aiming  vector  (-  a)  from  the  centroid  of  the 
area  of  coverage  is  not  too  extreme. 


2.  SYSTEM  GEOMETRY 

Radiation  from  a feed  (which  is  fixed  with  respect  to  the  ship)  passes  throu^  a hole 
in  a mirror  M and  falls  upon  a fixed  reflector  S,  which  reflects  it  back  toward  M.  The 
radiation  is  linearly  polarized,  and  S is  composed  of  wires  aligned  parallel  to  the  direction 
of  polarization.  The  mirror  acts  as  a half-wave  plate  and  rotates  the  angle  of  polarization 
by  90°,  so  that  the  radiation  passes  out  throu^  S after  reflection  by  M.  Let  a be  the  unit 
vector  defining  the  direction  of  the  radiation  incident  on  the  mirror.  Then  a is  fixed  to  the 
ship  and  points  into  the  ship,  and  its  negative  (-  a)  will  be  called  the  aiming  vector. 

Let 


p = unit  vector  normal  to  the  mirror,  pointing  outward; 


{ = “reflected”  unit  vector,  pointing  to  the  target. 


Then,  since  the  angle  of  incidence  on  the  mirror  equals  the 
an^e  of  reflection,  the  vectors  a,  p,  and  ( are  related  by 
(Fig.l) 


{ = a-2(a-p)p  (2.1) 


and 


Az. 

l{-a|  • 


a 


(2.2) 


Since  target  ranges  have  no  relevance  to  our  problem, 
we  can  identify  each  target  position  with  a point  on  a unit 
sphere,  whose  position  is  (.  The  corresponding  value  of  p 
(given  by  (2.2))  will  be  called  the  p-image  of  the  target. 

Let  be  an  arbitrary  collection  of  ^tar- 

get positions,  and  let  p^^\  ...,  p^^^  be  the  corresponding 
p-images.  The  steering  problem  is  to  find  the  closed  curve 
p = p(0  (0  < f < 7),  subject  to  certain  constraints,  which 


I 


Fig.  1 — Relation  of  the  vectors 
a,  p,  and  ( 
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joins  all  the  p-images  with  a minimum  update  time  T.  In  the  steering  problem  it  is  assumed 
that  the  value  of  the  aiming  vector  (-  a)  is  fixed  and  given.  The  aiming  problem  is  to  make 
an  intelligent  choice  of  a for  each  given  area  of  coverage. 

We  shall  fix  a ri^t-handed  rectangular  X1X2X3  coordinate  system  to  the  ship,  with 
corresponding  spherical  coordinates  r,  6,  and  0.  Ajs  Aown  in  Pig.  2,  the  mirror-scan  system 
is  attached  to  the  port  side  of  the  ship,  the  posi- 
tive Xi  axis  points  toward  the  bow,  the  positive 
axis  is  vertical,  and  the  coordinates  of  a point 
on  the  unit  sphere  (r  = 1)  are 


Xi=  and  cos 0 ' 
Xg  “ sin  0 sin  0 
X3  = cos  6 


(2.3) 


Since  a points  into  the  ship,  its  Xg  coordinate  is 
negative.  Similarly,  the  Xg  coordinates  of  p and  ( 
are  always  positive. 


3.  THE  CONTROL 

Mathematically,  the  mirror  motion  is  defined  by  a curve  p = p(0  lying  on  the  imit 
sphere,  that  is,  by  a description  of  the  motion  of  the  unit  normal  p.  Physically,  the  mirror 
is  moved  by  applying  torques  to  the  mirror  around  certain  instantaneous  axes  of  rotation, 
and  the  system  control  is  defined  by  specifying  the  total  torque  vector  u as  a function  of 
time:  u = u(t).  The  relationship  between  the  control  u and  motion  of  p is 

u»/pXp,  (3.1) 

ii^ere  / is  a certain  moment  of  inertia. 

Equation  (3.1)  is  a special  case  of  Euler’s  equation  for  the  motion  of  a rigid  body,  and 
in  its  derivation  two  assumptions  are  made: 

Al.  The  mirror  is  a disk,  a square,  or  some  other  body  having  certain  symmetry  prop- 
erties. For  a disk,  / is  the  moment  of  inertia  about  a diameter. 

A2.  The  torques  are  applied  in  such  a manner  that  the  mirror  motion  has  no  rotational 
component  around  the  axis  p,  because  the  existence  of  such  a component  would  imply  the 
expenditure  of  energy  to  produce  motions  of  the  mirror  which  have  no  effect  on  its  aim. 

It  turns  out  to  be  more  convenient  to  write  (3.1)  in  the  form 

u=/pX52p/5f^  (3.1*) 
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where,  by  definition, 

6^p/6t2  = p + |p|2p.  (3  2) 

The  derivation  of  (3.1*)  and  thus  of  (3.1)  is  given  in  Appendix  A.  The  mathematical  equiv- 
alence of  (3.1)  md  (3.1*)  is  a trivial  consequence  of  the  identity  p X p = 0,  which  allows 
any  scaler  multiple  of  p to  be  added  to  the  p term  in  the  right-hand  side  of  (3.1)  without 
changing  its  value.  The  replacement  of  p in  (3.1)  by  p + |p  l^p  and  thus  by  may 

therefore  appear  to  be  hij^ly  arbitrary,  but  it  is  justified  by  the  identity 

IpX  P 1=  IS^p/Sf^i,  (3.3) 

so  that  the  constraint  on  the  torque  magnitude  |u|  reduces  to  a constraint  on  |6^p/6f^  |.  To 
show  that  (3.3)  is  an  identity,  one  differentiates  the  identity  p • p = 1 twice  and  obtains 

p*p  = 0andp*p+  |p|2  = 0.  (3.4) 

Equation  (3.3)  is  a consequence  of  (3.4)  and  of  the  vector  identities 

a • (b  X c)  = - c • (b  X a) 


and 


a X (b  X c)  = b(a  ' c)  - c(a  • b). 
From  (3.2)  and  (3.4)  we  also  have 


|52p/5f2|2  = |p|2  _ |^|4  ^3  gj 

The  quantity  6^p/5f^  is  called  the  second  covariant  derivative  of  p = p(t),  and  in 
Appendix  B it  is  shown  that  S^pIBt^  is  the  tangential  component  of  p,  that  is,  the  projec- 
tion of  the  vector  p(f)  onto  the  plane  tangent  to  the  unit  sphere  at  the  point  p(f).  If  the 
curve  p ” p(f)  is  taken  to  be  the  description  of  the  motion  of  a particle  of  unit  mass,  then 
6 p/6f  is  equal  to  the  tangmtial  component  of  the  force  acting  on  the  particle.  In 
particular,  the  condition  = 0 is  precisely  the  condition  that  the  motion  be  inertial, 

that  is,  on  a great-circle  arc  with  constant  speed.  Hence  (3.1*)  has  the  following  conse- 
quence: if  the  applied  torque  u is  cut  off  at  a certain  time,  then  p = p(f)  will  thereafter 
move  along  a great-circle  route  with  constant  speed. 


4.  THE  STEERING  PROBLEM 

In  this  swtion  we  define  the  steering  problem  in  its  most  general  form.  The  con- 
straints of  this  general  problem  are  easy  to  express  and  hard  to  apply,  but  in  our  subsequent 
discussion  of  the  stop-go-stop  procedure  (section  5)  we  will  show  how  these  constraints 
reduce  to  a single  constraint  which  is  simply  and  directly  related  to  the  motion  of  p. 

The  first  constraint  can  be  stated  simply  as 

Iu|<a',  (4.1) 
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where  a'  is  the  largest  torque  that  can  be  applied  to  the  mirror.  This  constraint  is  related  to 
the  motion  of  p via  (3.1)  or  (3.1*)  and  can  be  written  in  one  of  the  following  equivalent 
forms,  each  of  which  will  hereafter  be  referred  to  as  the  constraint  (4.C1): 


or 


or 


Ip  X p|<a, 

|62p/6f2l<o, 

lp|2-|pl*<a. 


(4.C1) 


where  a = oi'll,  I being  the  moment  of  inertia  which  appears  in  (3.1)  and  (3.1*).  The  mathe- 
matical equivaloice  of  these  forms  is  a consequence  of  (3.3)  and  (3.5). 


The  second  constraint,  referring  to  the  beam  velocity  ^ (I)  as  the  beam  passes  throu^  a 
target,  will  be  expressed  as 


IS  (■>•2) 

where  0 is  the  largest  possible  beam  velocity  which  allows  the  required  number  of  radar 
pulse  hits  on  the  target.  Differentiating  (2.1)  to  express  (4.2)  in  terms  of  the  motion  of  p, 
we  get  the  constraint 


|(a  • p)p  + (a  • p)p  < 0/2.  (4.C2) 

Using  the  same  notation  introduced  in  section  2,  we  let  be  the  position 

vectors  of  N targets  and  let  p^^\  ...,  p^^^  be  the  corresponding  p-images.  We  let  { = {(I)  be  a 
curve  which  joins  the  points  ...» and  we  let  p = p(I)  be  the  p-image  of  this  curve, 
so  that  for  each  value  of  t,  p(f)  and  {(t)  are  related  by  (2.1)  and  (2.2).  We  can  now  state  the 
steering  problem  in  its  most  general  form: 

The  Steering  Problem.  Find  the  closed  curve  p = p(t)  which  satisfies  the  constraints  (4.C1) 
and  (4.C2)  and  joins  all  the  points  p(l), ...,  p^^^  with  minimum  update  time. 

Remark  4.1.  Later,  vhen  the  design  of  the  system  has  been  made  more  definite,  it  may  be 
necessary  to  replace  the  constraint  (4.1)  (and  hence  constraint  (4.C1))  by  a constraint 
imposing  separate  bounds  on  the  components  of  u about  certain  instantaneous  axes  of 
rotation.  A further  complication  mi^t  also  be  required  if  the  moments  of  inertia  about 
these  axes  shoxild  have  too  strong  a functional  dependence  on  position.  In  this  case  the 
constraint  (4.1)  would  be  only  an  approximation  but  a safe  one,  since  the  values  of  a'  and  a 
would  always  be  adjusted  so  as  to  guarantee  that  no  component  of  u exceeds  a certain 
amount,  and  the  estimation  of  the  update  time  T thus  obtained  would  be  greater  than  its 
true  value. 

Remark  4.2.  An  attractive  alternative  to  constraint  (4.C2)  would  be  the  “approximation” 

2|p|<0.  (4.C2A) 
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In  fact,  by  setting  a ' p = - cos  , it  is  easy  to  show  that  |a  * p | < sin  <//  Ip  |,  so  that  from 
(2.1)  one  obtains 

|{l2  - 4 |(a  • p)p  + (a  • p)p|2 
= 4[(a  • p)2  +(a  • p)^lpl2] 

<4[sin^  ^ +cos^  \I']1p1*" 

In  other  words 

l|l<2|p|.  (4.3) 

so  that  approximation  (4.C2A)  is  safe  in  tht  same  sense  as  in  the  previous  remark.  Moreover 
the  accuracy  of  the  approximation  l||  = 21p|  is  usually  quite  good;  it  depends  on  the  maxi- 
mum permissible  excursion  of  p from  - a and  therefore  depends  on  the  area  of  coverage.  For 
example,  if  the  area  of  coverage  is  a quarter-hemisphere,  then  |||  and  2 |p|  differ  by  less  than 
8 percent. 


5.  THE  STOP-GO-STOP  PROCEDURE 

As  before,  we  let  be  N target  positions  and  let  p^^^  ...,  p^^^  be  the  cor- 

responding p-images.  In  the  stop-go-stop  (SGS)  procedure  p = p(t)  moves  along  great-circle 
h arcs  in  a closed  path  joining  all  the  p^''^  and  comes  to  a complete  stop  at  each  of  the  p^-'X 

Hence  constraint  (4.C2)  is  automatically  satisfied,  and  we  will  now  show  how  constraint 
(4 .Cl)  reduces  to  a mathematically  simpler  form,  which  will  be  denoted  by  (5.CSGS). 

Our  first  task  is  to  determine  the  ciurve  p = p(f)  which  joins  two  points,  say  p^^^  and 
p^^\  in  the  shortest  time,  su^ect  to  constraint  (4.C1).  We  let  Z(l,  2)  be  the  great-circle 
distance  between  p^^^  and  p^^  (the  arc  length  of  the  short  great-circle  arc  joining  these 
points).  Hien 

Z(l,  2)  = cos'^  [p<^>  • p(2)] . (5.1) 

Let  X = x(s)  (0  < 8 < Z(l,  2))  be  the  short  great-circle  arc  joining  p^^^  to  p^^^  parame- 
terized by  s = arc  length.  Then  x = x(8)  is  (geodesic)  motion  along  a great-circle  arc  with 
constant  unit  speed,  and  the  acceleration  vector  is  a unit  vector  which  is  always  directed 
toward  the  center.  Hence,  letting  primes  denote  differentiation  with  respect  to  s,  (')  = 

(d/ds),  we  have 


Ix'l  = 1,  (5.1) 

x'  • x"  = 0.  (5.21 
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(We  recall  from  (3.2)  that  by  definition  8^x/8s^  = x"  + Ix'Px,  so  that  8^x/88^  = 0 by 
virtue  of  (6.1)  and  (5.3);  that  is,  the  motion  x = x(s)  w geodesic.) 

Now  let 


P(f)  = x[s(0]  (5.4) 

be  a general  motion  along  the  great-circle  arc  whose  velocity  |p|  = Idp/dtl  is  not  necessarily 
constant.  Differentiating  (5.4),  we  get 


p = sx'  (5.5) 

and 

p = 8x'  + 8^x".  (5.6) 

But  p(f ) is  merely  a reparameterization  of  x(s),  so  that  from  (5.3)  we  have  x"  = - p.  Also, 
from  (5.2),  x'  is  a unit  vector  perpendicular  to  p.  Referring  to  (5.6),  we  therefore  see  that 

8x'  = tangential  component  of  p = 6^p/6f^, 


so  that 


Ip  X pi  = Isp  X x’l  = Is  |. 
Hence  constraint  (4. Cl)  reduces  to 


l«l<a,  (5.CSGS) 

where  a is  the  same  constant  as  appears  in  (4. Cl). 

The  problem,  of  finding  the  curve  p = p(f)  joining  p<^)  to  p^^^  in  minimum  time  T sub- 
ject to  constraint  (4.C1)  as  well  as  the  conditions  p(0)  = p(r)  = 0,  has  now  been  reduced  to 
the  problem  of  finding  the  piecewise  smooth  function  s = s(0  for  which 

(5.CSGS)  is  satisfied  on  each  smooth  piece, 

8(0)  = 0 and  8{T)  = Z(1,2), 

8(0)=8(D  = 0, 


and 


T is  minimum. 

This  is  a standard  problem  in  control  theory  and  has  a well-known  solution:  one 
applies  a maximum  acceleration  a for  half  the  time  and  a maximum  deceleration  - a for  the 
other  half.*  Hence  we  have 


•E.B.  Lee  and  L.  Markus,  Foundations  of  Optimal  Control  Theory,  Wiley,  New  York,  1967,  Ch.  1. 
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or,  solving  for  T ^ 


Tmin  = 2) la,  (5.7) 

TTiere  are  {N  - l)!/2  SGS  paths  joining  the  N points  and  we  now  consider 

tht^  problem  of  finding  the  one  which  is  optimal. 

Consider  for  example  the  configuration  of  four  points  shown  schematically  in  Fig.  3. 
The  points  are  labeled  1,  2,  3,  and  4,  and  there  are  three  SGS  paths  labeled  1234, 1324,  and 
1342.  For  example,  the  path  1234  is  the  path  which  takes  point  1 to  2,  2 to  3, 3 to  4,  and 
4 to  1.  These  three  paths  are  the  only  three  essentially  different  SGS  paths  in  this  case, 
since  all  the  other  paths  differ  from  these  only  in  their  starting  point  or  orientation.  Thus, 
for  our  purposes,  the  path  labeled  1234  is  the  same  as  paths  labeled  2341, 1432,  etc. 

More  generally,  for  a configuration  of  N points  we  can  always  label  the  starting  point 
as  1,  so  that  theiV-tuple  (1,  (2«  i'si  >7/)  throu^  all  the  different  path  levels  as  the 
(N  - l)-tuple  (i2f  13.  •••.  nms  throu^  all  the  permutations  of  the  {N  - 1)  integers  {2,  3, 
...,  N} . In  fact,  each  path  is  represented  twice,  once  for  each  orientation,  so  that  there  are 
only  (N  -1)1/2  essentially  different  SGS  paths  joining  N points. 

We  let  Z(i,j)  be  the  arc  length  of  the  great-circle  arc  joining  p^'Ho  ^ and  let  T{i,j) 
be  the  transit  time  tor  the  motion  of  p(t)  from  to  p^^^  the  motion  being  “optimal”  as 
previously  described  in  this  section.  Then  as  in  (5.1)  and  (5.7)  we  have 

Z(i,  j)  = cos-1  [p(i)  . p(y)]  (5  3) 

and 

T{iJ)-2y/2{i,i)la,  (5.9) 

and  the  update  time  T corresponding  to  the  path  (1,  i2>  >3  > —t  ) is 

T=  + V^«‘2.«‘3)  + - + V'ZO'n-i.'n)  + s/ZHnA)).  (5.10) 

V« 


Fig.  3 — Paths  joining  four  points 
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Hence  the  optimal  SGS  path  is  the  one  which  minimizes  the  ri^t-hand  side  of  (5.10).  In 
brief,  we  shall  say  that  the  optimal  SGS  path  minimizes  Hy/Z. 

For  the  moderately  small  number  of  targets  considered  in  this  report  (JV  ^ 7)  the  opti- 
mal SGS  path  can  be  found  by  the  straightforward  enumeration  of  all  - l)!/2  cases.  The 
general  problem  of  minimizing  a “cost  function”  such  as  (5.10)  over  a class  of  paths  is 
known  as  the  traveling-salesman  problem,  and  algorithms  have  been  constructed  which 
make  the  search  for  the  optimal  path  more  efficient.  However  these  algorithms  are  all 
inefficient  in  the  sense  that  the  computational  time  varies  exponentially  with  N.*  Much 
attention  has  been  given  to  the  search  for  optimal  and  near-optimal  solutions  to  the  prob- 
lem of  minimizing  2Z.  However,  simple  examples  can  be  constructed  which  show  that 
paths  which  minimize  £Z  do  not  necessarily  minimize  "Zy/Z,  and  conversely. 


6.  THE  AIMING  PROBLEM 

Figure  4 shows  a half -hemisphere  with  four  “extreme”  points  on  the  boundary  and 
their  coordinates.  The  coordinate  system  shown  here  is  the  same  as  that  shown  in  Fig.  2. 

As  5 varies  over  the  entire  half -hemisphere,  its  p-image  p will  vary  over  a smaller  region,  and 
the  range  of  p depends  on  the  value  of  the  aiming  vector  a.  This  is  shown  in  Fig.  5 for  three 
values  of  a.  The  outer  curve  is  the  boundary  of  the  half-hemisphere,  as  seen  by  a distant  eye 
having  spherical  coordinates  (Fig.  2)6  = 60°  and  <t>  = 90°.  The  inner  curve  is  the  boundary 
of  the  range  of  p,  and  it  is  marked  off  into  two  subregions  I and  II  which  are  the  p-images 
of  the  quarter-hemispheres  whose  extreme  points  are  (1,  2, 4)  and  (2, 3,  4)  respectively. 

Table  1 shows  how  the  great-circle  distances  between  the  p-images  of  the  extreme 
points  varies  with  three  values  of  the  aiming  vector  a.  TOe  arc  lengths  Z(i,  /),  in  degrees,  are 
computed  according  to  (5.8),  in  which  the  values  of  p^*^  are  computed  according  to  (2.2). 
Z(l,  3)  is  always  90°,  because  the  position  vectors  J of  points  1 and  3 are  the  negatives  of 
each  other,  so  that,  from  (2.2),  p^^^  • p^^^  = 0 for  every  value  of  a.)  The  results  shown  in 
this  table  and  other  calculations  suggest  that  the  second  value  of  a listed  is  optimal  when 
the  area  of  coverage  is  the  entire  half-hemisphere  and  that  the  third  value  is  optimal  when 
the  area  of  coverage  is  restricted  to  the  quarter-hemisphere  (1,  2,  4). 


X3 


Fig.  4 — Coordinates  of  four  "extreme” 
points  on  a half-hemisphere 


•H.R.  Lewis  and  C.H.  Papadimitriou,  "The  Efficiency  of  Algorithms,”  Scicnfi/ic  American  238  (No.  1), 
96-109  (Jan.  1978). 
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Fig.  5 — Ranges  of  p-images  p for  vajrioua  aiming  vectora 
divided  into  regiona  1 and  II  (the  ends  of  the  subdividing  ai 
and  4) 


-a  « (-1/73,1/^,  1/73)  j 

I,  with  the  ranges  of  p sub-  i 

s being  p-images  of  points  2 ! 


7.  CONCLUSIONS 

7J.  Calculations  of  Optimal  SGS  Paths  and  Update  Times 


Optimal  SGS  paths  and  update  times  were  computed  in  the  manner  described  in 
section  5,  with  the  current  design  value  of  the  constant  a (defining  the  constraint)  given  by 

a = (lOOv'5)  radians/s^  « 8103  deg/s^. 

I 

[ 

Optimal  paths  and  update  times  T were  computed  for  various  configurations  of  up  to 
seven  points  and  various  values  of  a.  Some  examples  of  our  calculations  are  given  in  Fig.  6 
for  the  three  values  of  a listed  in  Table  1.  The  half-hemisphere  shown  in  Fig.  4 is  shown  in 
[ Fig.  6 as  seen  by  an  observer  who  is  in  the  X2X^  plane  at  large  and  positive  and  slightly 

[ negative  X2  • The  target  positions  ( are  indicated  by  heavy  dots,  and  N is  the  number  of 

targets. 


7 Summary  of  Observati<ms  Concerning  Calculated  Results 


t 


The  results  of  these  calculations  can  be  summarized  in  the  following  observations : 

i 


Observation  1.  The  Mirror  Scan  system  can  simultaneously  track  at  least  seven  targets 
in  a half-hemisphere  with  an  update  time  of  less  than  1 second. 

Observation  2.  The  value  of  the  aiming  vector  - a has  only  a sli^t  effect  on  the  effici- 
ency of  the  system;  for  the  three  values  of  a listed  in  Table  1,  the  variation  in  update 
time  was  about  5%  or  less. 


Table  1 — Variation  of  Z(i,  J)  with  a 


- a 

X(l.  2) 

Z(l,3) 

^(1.4) 

1 

Z(2.3) 

^(2,4) 

Z(3,4) 

[0.1,0] 

60.0° 

90.0° 

60.0° 

45.0° 

[0,  lA/J,  l/v^2] 

49.2° 

90.0° 

49.2° 

49.2° 

49.2° 

[-l/v^.l/v/J,lV3] 

46.9° 

90.0° 

46.9° 

52.2° 

52.2° 
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Appendix  A 

DERIVATION  OF  EQUATION  (Sa) 

As  explained  in  section  3,  equations  (3.1)  and  (3.1*)  are  mathematically  equivalent, 
and,  to  prove  (3.1*),  we  shall  first  express  p in  terms  of  certain  angular  accelerations  and 
velocities. 

As  shown  in  Fig.  A,  let  and  62  be  two  unit  vectors  which  are  fixed  to  the  mirror  M 

in  such  a way  that  {e^ , 62 , p}  form  a ri^t-handed  triple: 


ej  X 62  = P,e2  X p = ei,p  X Bi  =02.  (Al) 

Since  61  ’ 01  = 1,  we  have  0i  * ^i  = 0,  so  that  0i  has  no  component  in  the  01  direction. 
Also,  from  assumption  A2  in  section  3,  it  follovrs  that  Oi  has  no  component  in  the  02 
direction.  Similar  remarks  hold  for  02,  so  that  one  can  define  angular  velocities  Ui  and 
satisfying 

01  * - WiP  and  62  “ - W2P* 

From  equations  (Al)  and  (A2)  we  get 


p • WiBi  + <I>2«2 


and 


p-cii0i  + CO202  - (wf  + 

Since  8^p/5t^  is  by  definition  the  tangential  component  of  p,  it  follows  that 


(A3) 

(A4) 


6*p/5f*  - Wj0i  + W2«2‘ 


(A6) 


P 


Fig.  A — Vectors  relative  to  a disk- 
shaped  mirror 
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We  will  now  compute  the  torque  u.  It  will  be  convenient  to  think  of  the  mirror  as 
being  composed  of  particles  all  having  the  same  unit  mass.  Let  r,-  be  the  position  vector  of 
the  ith  particle,  and  set 

I'*  “ ^ 

Then,  by  definition, 

u ■ Af  X fj.  (A6) 

Each  particle  has  “local"  coordinates  x}  and  xf  defined  by 

r,.  = xjBi  + xfeg, 

and  xj  and  xf  are  constant,  since  the  moving  axes  and  82  are  fixed  to  the  mirror.  From 
(A6)  one  ultimately  gets 

u = - [/i2^  ^22^1®!  ^12^1®2 

[^12(^1  “ ^2^  " (^11  ~^22)^1^2lP» 

where 

^12  E ^22  = E (*?)*• 

i I i 

Rrom  the  symmetry  properties  of  the  disk  we  have 

Ii2  ~ 0 and  /ji  = /22  °°  f» 

where  / is  the  moment  of  inertia  about  a diameter.  Hence  for  the  disk  we  have 

u * /(-  ^6]^  + ^62) 

-/(«2P  X 62  +^p  X 81), 

and  firom  (A5)  we  get  (3.1*): 

u = /p  X 6*p/6t^. 

ii 

( 


! 


Appendix  B 

COVARIANT  DERIVATIVES 


.i 


[' 

t 


I 


► 


Let  p * p(t)  be  a curve  on  the  unit  si^eie  S.  We  want  to  diow  that  the  right4uuid  side 
of  (3.2)  is  the  tangential  component  of  p.  Since  the  vector  p is  normal  to  S at  the  point  p, 
we  have 

normal  component  of  p ~ (p  * p)p. 

Differentiating  the  identity  p * p = 1,  we  get 

p • p = 0 and  p • p + IpI*  =*  0. 


Hence 


normal  component  of  p - Ipl^p, 


and  it  follows  that 


tangential  component  ~ p - normal  componoit 


* P + Ipl^p. 


The  preceding  formulas  apply  only  to  spheres.  More  generally  let  S be  any  surface  in 
Euclidean  three-space  defined  by  the  parametric  equations 

X'“X»(yl,y2),  i-1.2.3. 

where  and  y2  local  surface  coordinates.  Let  p(t)  = X(y(0)  be  a smooth  curve  on  S. 
Then  the  second  covariant  derivative  of  p(t)>  , is  again  defined  as  the  tangential 

component  of  p(t)  and  is  given  by 


where 


8^ 

at* 


ax  ^ 6^ 
3yi  at* 


ax 

ay*  ’ 


a*y° 

at* 


. yO  + 


L 

P.y 


and  where  are  the  so-called  Christoffel  symbols.  Details  can  be  found  in  almost  any 
elementary  text  in  differential  geometry  (for  example,  A.  J.  McConnd,  Applications  of 
Tensor  Analysis,  Dover,  New  York,  1967). 
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